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We have performed a photoemission study of Ti40 7 around its two transition temperatures so as to 
cover the metallic, high-temperature insulating (bipolaron-liquid), and low-temperature insulating 
(bipolaron-crystal) phases. While the spectra of the low-temperature insulating phase show a finite 
gap at the Fermi level, the spectra of the high-temperature insulating phase are gapless, which 
is interpreted as a soft Coulomb gap due to dynamical disorder. We suggest that the spectra of 
the high-temperature disordered phase of Fe304, which exhibits a charge order-disorder transition 
(Verwey transition), can be interpreted in terms of a Coulomb gap. 

PACS numbers: 71.23.-k, 71.30.+h, 72.80.Ga, 79.60.-i 



Since Mott jjj proposed the idea of variable-range 
hopping and minimum metallic conductivity for disor- 
dered systems and Anderson 0| raised the concept of 
localization due to disorder, physical properties of dis- 
ordered solids have been extensively studied. Influence 
of Coulomb interaction on the electronic density of states 
(DOS) near the Fermi level (Ep) of disordered systems is 
one of the most fundamental issues to be clarified. Efros 
and Shklovskii || proposed that in disordered insulators 
long-range Coulomb interaction opens a soft Coulomb 
gap (SCG), whose DOS is proportional to (E — E F ) 2 . 
So far, there have been tunneling spectroscopic confir- 
mations of the SCG in some disordered systems such as 
doped semiconductors [Q. Davies and Franz || pointed 
out the possibility of an SCG opening in the photoe- 
mission spectra of Na^Ta^Wi-yOs but the experi- 
ments did not have sufficient energy resolution to criti- 
cally address this problem. Another important issue is 
how short-range order (SRO) affects DOS near Ep, that 
is, how the electronic structure evolves when a charge 
ordering develops from disorder to SRO to long-range 
order. 

Ti407 is a suitable system to study the above prob- 
lems: It undergoes successive phase transitions with de- 
creasing temperature from a metal to a charge-ordered 
insulator via an insulator with SRO |Q. It is a system 
with nominally 0.5 3<i electron per Ti, allowing two pos- 
sible valence states of Ti 3 + (3d 1 ) and Ti 4+ (3d ), and at- 
tracted particular attention in 1970's as a system where 
bipolarons, or singlet pairs of two polarons, are formed 
in real space @-|jJ. Above Tmi — 154 K, the system 
is in the metallic (M) phase and the Ti valence is be- 
lieved to be uniform 3.5+ as shown in Fig. [l] (a). That 
is, the electrical resistivity p(T) is metallic with Pauli- 
paramagnetic x(T)- With decreasing temperature, sin- 



glet Ti 3+ -Ti 3+ pairs, namely bipolarons, are formed, re- 
sulting in the metal-to-insulator transition at Tmi with a 
steep increase in p(T) by three orders of magnitude. At 
the same time, x(T) almost vanishes, reflecting the for- 
mation of the singlet bipolarons. We refer to this phase as 
the high-temperature insulating (HI) phase. Because the 
bipolarons are dynamically disordered in this phase as 
has been established by EPR studies Jl(|, the HI phase 
may be called a bipolaron liquid. Subsequently, bipo- 
larons become ordered below Tn ~ 140 K as shown in 
Fig. Q (a), with a further rise in p(T) by three orders of 
magnitude while x(T) remains unaffected. We refer to 
this phase as the low-temperature insulating (LI) phase. 
The transition from the liquid to the crystal of bipolarons 
is a kind of Verwey transition pljji^] and the bipolaron- 
liquid formation in the HI phase is a kind of a SRO of 
charge carriers. The Verwey transition Jl3| was originally 
found for Fe 3 4 at T v ~ 120 K. Analogous to Ti 4 7 , 
Fe30 4 has two possible valence states of Fe 2+ and Fe 3+ 
which are ordered below Ty and are disordered above 
Ty, leading to the characteristic jump in the electrical 
resistivity. 

In this Letter, we have studied the single-particle exci- 
tation spectra of the three phases of Ti 4 7 by means of 
photoemission spectroscopy (PES) with high energy res- 
olution. To add to the variation of the electronic struc- 
ture of Ti 4 7 as a function of temperature across Tmi 
and Trj, we focus on general aspects of the single-particle 
excitation spectra of dynamically disordered systems and 
then made comparison with the PES result of Fe3 4 . 

Single crystals of Ti 4 7 were grown by the vapor trans- 
port method [l^j . PES measurements were performed 
using an OMICRON hemi-spherical analyzer and a He 
lamp (He I: hv —21.2 eV). The energy calibration and 
the estimation of the instrumental resolution were done 
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by measuring the Fermi edge of Au evaporated on the 
sample. The energy resolution was set ~50 meV. Sam- 
ples were cleaved in situ. This always gave an irregular 
surface, an assembly of randomly orientated small facets. 
As the analyzer had an acceptance angle of ±8°, the ob- 
tained spectra can be regarded as angle-integrated PES 
spectra. The measurement temperature was controlled 
within the accuracy of ±0.2 K. The base pressure in the 
spectrometer was less than 1 x 10 -10 Torr. Below, we 
show results reproducible for several cleaves and obtained 
within a few hours after cleavage. 

The position and the width of the O 2p band (not 
shown) show almost no temperature dependence across 
the two phase transitions as in the previous photoemis- 
sion study jl(|. In contrast, the PES spectra in the Ti 
3d band region show strong temperature dependence as 
shown in Fig. |l| (b), being consistent with the previous 
results for the M and LI phases with a lower energy res- 
olution [fl6|| . Here, the sample was first cooled and then 
heated as indicated in the figure, so as to cross twice each 
transition temperature. Two observations are remarked. 
First, as superimposed in the bottom panel, there are 
three kinds of spectra, which reflect the electronic struc- 
ture of the M, HI, and LI phases of Ti4C>7 as discussed 
below. Second, the 142 K spectra show different spectral 
lineshapes between cooling and heating due to thermal 
hysteresis between the HI and LI phases, as has been ob- 
served in the electrical resistivity, thermopower [Q, and 
EPR measurements [jul. In Fig. |l| (c), we have plotted 
the temperature dependence of the integrated PES in- 
tensity within 0.5 eV of Ep normalized to the intensity 
integrated from Ep to binding energy Eb = 1 eV. This 
hysteretic behavior was reproducible for several thermal 
cyclings across Tmi and Tjj. 

Let us discuss the spectra of each phase in more de- 
tail. As shown in Fig. [I] (d), M-phase spectra show a 
weak but finite Fermi edge, reflecting its metallic na- 
ture. The spectrum shows a broad maximum centered at 
Eb = 0.75 eV, around which most of spectral weight is 
distributed. A similar broad feature has been observed 
at Eb ~ 1-1.5 eV in the spectra of three-dimensional 
titanium and vanadium oxides with 3d 1 configuration, 
LaTiC>3 and YTiC>3 jl7| , and has been interpreted as the 
incoherent part of the spectral function accompanying 
the coherent quasi-particle excitations around Ep. The 
coherent part in Ti4C>7 is hardly separable from the inco- 
herent part because of their strong overlap, leading to the 
pseudo-gap-like behavior at Ep. This implies that the M 
phase of Ti4C>7 is a strongly correlated metal, where the 
motion of conduction electrons is largely incoherent, re- 
sulting in the weak coherent part. In going from the M 
to the LI phase, the incoherent peak at Eb = 0.75 eV be- 
comes sharper and a clear gap of order ~0.1 eV is opened 
as seen m Fig. (d). The gap opening is attributed to 
the bipolaron ordering in this phase. The Hl-phase spec- 
tra take an intermediate lineshape between the LI and 
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FIG. 1. (a) Schematic crystal structures of Ti407 which 
show only the chains parallel to the pseudo-rutile c-axis trun- 
cated by shear planes indicated by thick lines [B. In the LI 
phase, the bipolarons are ordered. In the M phase, all the 
Ti sites have a uniform valence of 3.5+. (b) PES spectra of 
Ti407 taken by varying the temperature from the bottom to 
the top. (c) Normalized PES intensity near Ef (see text). 
The solid line is guide to the eye. (d) The M-, HI-, and 
Ll-phase spectra around Ef plotted on an expanded scale. 



M phases: The 0.75 eV peak is somewhat broader than 
in the LI phase. The spectra show neither a Fermi edge 
nor a real gap. Indeed, the PES intensity vanishes only 
at Eb = 0.0 eV as shown in Fig. [I] (d). Figure || (a) 
shows the same spectra plotted against E B sgn(Es). The 
Fermi edge in the M-phase spectra is now clearer while 
the LI spectra, which show a "hard gap" , are concave 
around Ep. Most remarkably, the Hl-phase spectra form 
an almost straight line from E B = 0.00 to ~ 0.08 eV 2 
(0.0 < Eb ^ 0.3 eV), meaning that the spectra show a 
power-law behavior with the exponent of ~ 2. 

To be more quantitative, we have performed a line- 
shape analysis for the LI- and Hl-phase spectra tak- 
ing the experimental resolution into account. The line- 
shape was assumed to be of the form I(Eb,T) = 
I (EB,T)+I hg (E B ), where I (E B ,T) is the intrinsic part 
expressed by I (E B ,T) = a 2 E 2 B f(E B ,T). f(E B ,T) = 
[exp(— .Eb /k B T) + 1] _1 is the Fermi-Dirac distribution 
function at T although the finite-temperature effect due 
to f(EB,T) is negligibly small because I o (0,T) = 0. 
Ibg(EB) = bo represents a constant background of the 
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FIG. 2. (a) PES spectra of Ti 4 7 plotted against 
E% sgn(_Es) with the fitted curves. The experimental spec- 
trum and the fitted curve plotted on the log-log scale for the 
HI phase (b) and the LI phase (c). 



PES spectra. For the LI phase, in order to represent the 
finite gap, we assumed I (E B ,T) — a 2 (E B ~E ) 2 Q(E B — 
Eq), where Q{x) is a step function and Eq was allowed to 
take a finite positive value. The instrumental resolution 
was included through the convolution of I(E B ,T) with 
a Gaussian of FWHM 50 meV. a 2 and bo were treated 
as adjustable parameters. As shown in Fig. || (a), the 
result of the fitting is satisfactory, especially for the HI- 
phase spectra for —0.1 < E B < 0.3 eV. One Hl-phase 
spectrum is shown in Fig. || (b) on the log-log scale to 
emphasize the E g dependence and the excellent fit for 
small E B . This makes a clear contrast to the Ll-phase 
spectrum shown in Fig. |^ (c), which shows a finite gap 
of£ . 

We propose that the above behavior of the HI phase 
can be explained as an SCG in the DOS, N{E B ) — 
^(■*k) 3 E B where k is the dielectric constant. For 
Ti 4 7 , the magnitude of the SCG A c = e 3 (A /K 3 ) 1/2 , 
where No is the noninteracting DOS at Ep, is estimated 
to be -0.2 eV if we take N - 0.01 eV" 1 !" 3 from the 
free-electron-like DOS around Ep and k — 10. As the 
estimated Ac and the observed soft gap have the same 
order of magnitude, we may conclude that the Hl-phase 
spectra are consistent with the opening of an SCG. Here, 
it should be noted that, even if there exists substantial 
SRO, i. e., bipolaron- liquid formation, in the HI phase 
of Ti407 the system is sufficiently disordered over long 
distance for an SCG to appear around Ep. That is, the 
length scale of the SRO, which is equal to the Ti-Ti dis- 
tance —3 A 0, is sufficiently shorter than the critical dis- 
tance Rc — e 2 /nA c ~ 8 A for the SCG to survive. The 
SRO would be reflected on high-energy spectral features, 
e. g., the sharpening of the E B — 0.75 eV peak in the 
HI phase. Experimentally, an SCG was observed in tun- 
neling spectroscopy measurements of the doped semicon- 



ductor Si:B Q), whose gap was as small as 0.75 meV due 
to low Nq. As for PES measurements, besides the afore- 
mentioned Na^TayWi-yOa, possible existence of SCG in 
Ti407 and Fe304 was suggested ||, while no quantitative 
analysis on the experimental spectra has been performed 
so far. In this sense, the present result I oc E 2 B for the 
HI phase of T14O7 is the first clear indication of an SCG 
using PES. Since our measurements were performed on 
cleaved surfaces, there are no extrinsic disorder effects 
induced by scraping. 

It is worth comparing T14O7 in the HI phase with 
Fes04 above Ty — 120 K because of the analogy be- 
tween the two materials pointed out repeatedly |p.2| . For 
this purpose, we reanalyzed the PES spectra of Fe304 
reported by Chainani et al. |]l4| in the context of an 
SCG. Figure || (a) shows the PES spectra of Fes04 mea- 
sured with the resolution of ~70 meV. They claimed 
that the intensity at Ep in the metallic phase evolves 
as the temperature increases. In Fig. || (b), we have 
replotted the spectra against E B sgn(E B ), which makes 
their argument clearer. More interestingly, we find that 
the 140 K (> TV) spectrum almost falls on a straight 
line at E B > 0.07 eV. In contrast, the 100 K spec- 
trum is slightly concave, signaling the opening of a hard 
gap. Fe304 is, however, different from T44O7 in that 
the 140 K spectrum, which would correspond to the 
Hl-phase spectra of T14O7, does not vanish at Ep but 
shows a finite Fermi edge. We have modeled the spectra 
I(E B ,T) above T v as I(E B ,T) = I {E B ,T) + I bg {E B ), 
where I (E B ,T) = {a 2 E 2 B + a )f{E B ,T) (a : finite 
DOS at Ep). The 100 K spectrum was assumed to be 
I (E B ,T) = a 2 (E B -E ) 2 9(E B - E ) as in the LI phase 
of Ti 4 07. Ib g (E B ) = b\E B + b represents the linear 
background of the PES spectra |l8[ ]. The instrumental 
resolution was also included. 

The result of the fitting was satisfying between E B = 
—0.15 eV and 0.15 eV as shown by solid curves in 
Figs. I (a) and (b). The dashed curves m Fig. I (a) 
for T > 140 K represent the intrinsic DOS cut-off at 
Ep {a 2 E 2 B + ao)f(E B ,0). The most remarkable point 
is that the finite intensity at Ep evolves systematically 
with increasing temperature. Figure]^ (c) shows the ratio 
ao/a 2 between the intensity at Ep, which contributes to 
the "metallic" conductivity, and the magnitude of the E B 
term, which represents the dynamical disorder. The elec- 
trical dc conductivity data of Fe304 [[l9| is superimposed 
in the figure. The ao/a 2 values at 140, 200, and 300 K 
form a straight line which extrapolates to at 119±5 K 
~ TV. This observation, namely, ao/a 2 oc T—Ty strongly 
indicates that in oc E B just above Ty. As the tempera- 
ture goes up above Ty, the finite intensity at Ep grows up 
in proportion to T—Ty. Park et al. [M reported that the 
spectra of F83O4 at 130 K show no Fermi edge, which is 
naturally understood as due to the small ao/a 2 . We pro- 
pose that an SCG exists in F83O4 just above Ty and con- 
tinuously evolves into a pseudogap well above Ty , reflect- 
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FIG. 3. (a) PES spectra of Fe 3 4 (hv = 21.2 eV) 
taken from Ref. Solid, dashed, and dotted-dashed 

curves represent the fitted curves, (aiE% + ao)f(EB,0), and 
<12(Eb — Eo) 2 8(Eb — Eo), respectively, (b) The same spectra 
plotted against E% sgn(Z?e). (c) Ratio an/a2, which tends to 
vanish at T ~ Ty (see text). The electrical conductivity of 
Fe304 is also shown by a dashed line 

ing the crossover from the semiconducting (dp/dT < 0) 
behavior just above Ty to the metallic (dp/dT > 0) 
one well above Ty. Just above Ty, charges are disor- 
dered with SRO, namely, the system is in a Wigner glass 
state Jjj as is the HI phase of T14O7. Here, it should 
be noted that the M-phase spectra of Ti 4 07 near Ep 
(E B < 0.2 eV) may be fitted to (a 2 E B + a )f(E B ,T) 
like the spectra of Fe3C>4 above Ty. However, a /a,2 = 2- 
4xl0~ 2 cV~ 2 is much larger than that in Fe3C>4 at 300 
K, reflecting the strongly first-order transition between 
the HI and M phases in T14O7. 

In conclusion, we have performed a PES study of 
TI4O7 covering its LI, HI, and M phases. The spectra 
of the Ti 3d band show peculiar temperature-dependent 
spectra characteristic of the three phases, among which 
the Hl-phase spectra are gapless and can be fitted to Tvfj 
near Ep. We interpret this as an SCG resulting from 
disordered bipolarons. By reanalyzing the PES spectra 
of Fe3C>4, an SCG was also found just above Ty, in- 
dicating the significant role of disorder and long-range 
Coulomb interaction. With increasing temperature, the 
SCG is found to continuously evolve into a pseudogap 
as the metallicity is gradually established, whereas the 
first-order HI-to-M transition in T14O7 precludes the ob- 
servation of the corresponding continuous change. Pre- 
sumably, similar SCG behavior may be observed in other 



systems such as P-Na x Y 2 05 M. The very recent scan- 
ning tunneling spectra of hole-doped manganites above 
T c [Ml may be interpreted in the same way as Fe 3 04 
above Ty. 
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